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Abstract 

A new QCD analysis of Deep Inelastic Scattering (DIS) data is presented. All 
available neutrino and anti-neutrino cross sections are reanalysed and included in 
the fit, along with charged- lepton DIS and Drell-Yan data. A massive factorisa- 
tion scheme is used to describe the charm component of the structure functions. 
Next-to- leading order parton distribution functions are provided. In particular, 
the strange sea density is determined with a higher accuracy with respect to other 
global fits. 



1 Introduction 



In Deep Inelastic Scattering (DIS) processes a neutral or charged lepton i inter- 
acts with a nucleon N yielding a lepton i' and a set of undetected hadrons X 
in the final state. The kinematics of this process, i{k) + N{p) £'{k') + X, 
is determined by two independent variables, besides the energy of the incoming 
lepton. One usually chooses them among the four Lorentz invariants 

Q'^-q' = -{k-kr,x = ^,y = ^, W' = {q + p)\ 

2p ■ q p ■ k 

which are obtained experimentally by measuring the momentum, the direction of 
the scattered lepton and the initial momenta k, p. 

Studies of DIS processes, both on the experimental ^ and on the theoretical 
side 0, have shed light on the nucleon internal dynamics. In particular, powerful 
tests of perturbative Quantum Chromo Dynamics (pQCD) have been carried 
out, consisting of global analyses performed on a large class of DIS observables. 
Among these the structure functions measured in neutrino (anti-neutrino) DIS 
play a major role for the determination of the fiavour distributions, in particular 
the valence and the strange sea densities. 

What experiments directly measure is a differential cross-section from which 
the structure functions are extracted by means of a theoretical analysis. This 
includes the application of electroweak radiative corrections, the determination 
of R = ai/ax, the incorporation of possible nuclear effects, etc ... In particular, 
the analysis of neutrino data is a difficult and delicate procedure (see for a 
clear and detailed description), and so far only a small part of the information 
accumulated in z/(z/)DIS has been exploited in the QCD parametrisations. 

The purpose of this paper is to present a new global analysis of DIS data 
which includes the available u and u cross section measurements^, besides the 
structure function data collected in charged-lepton DIS experiments. We resort 
directly to the i>{u) DIS differential cross sections, avoiding to use the neutrino 
structure function results. The latter are the product of a preanalysis which may 
be (and often indeed is) based on theoretical assumptions different from those 
of the global fit that one is performing. Thus, for a full consistency, we use 
only cross section data. This limits our neutrino (and anti-neutrino) data set to 
the BEBC (Hydrogen target) [^, CDHS (Hydrogen and Deuterium targets) 
and CDHSW (Iron target) § measurements. The CHARM § and CCFR § 
experiments do not provide cross sections but only structure functions and hence 
their results are not included in our analysis. 

The problem with the BEBC, CDHS and CDHSW data is that they cannot be 
used in the form they were published, since the electroweak radiative corrections 
were either incomplete or not applied at all and/or the bin centre corrections were 

-'^This type of analysis was proposed long ago by M.W. Krasny 
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not performed. Thus we had to reevaluate the neutrino cross sections to take all 
these corrections into account. This is the preliminary step of our analysis. Since 
most of the v (z/) DIS data come from nuclear targets, nuclear corrections must 
also be applied. 

Besides the neutrino data, the structure function measurements from charged- 



lepton DIS experiments (NMC |Tg, BCDMS [|TT], HI [g]) are included in our 
fits. The only non DIS data we use are the Drell-Yan [l^, |l5l data which 
constrain the light sea. 

Our fitting procedure is designed in such a way to take properly into account 
the experimental uncertainties and the correlations among them, which are known 
to affect the slopes of the structure functions and ultimately the determination 
of the parton densities . 



An important feature of the QCD analysis presented in this paper is the ac- 
curate treatment of the charm contribution to the structure functions. A massive 



factorisation scheme is used, the so-called Fixed Flavour Scheme . It is known 
[0 that for a precise extraction of the strange sea density, charm mass effects 
cannot be neglected and have to be correctly incorporated. 

In extracting the parton distributions, the neutrino data (in particular the 
high-statistics CDHSW data) add a great quantity of information to that com- 
ing from charged-lepton DIS. The latter is insufficient to constrain all flavour 
distributions, being essentially limited to one observable, F2. Charged-current 
DIS provides four more independent combinations of parton densities, Fg^, Fg^, 
xFg^, xF!^. As a consequence, the abundance of neutrino data in our fit ensures 
an excellent accuracy in the determination of the flavour densities. 

A special emphasis will be given to the strange sea density. Due to the lack 



of data able to constrain it, in the existing fits |T9|, ^ this distribution is tightly 
related to the non-strange sea distributions and essentially borrowed^ from the 
CCFR extraction |^2[. Clearly, this is not a consistent procedure. Here we present 
the first fully consistent determination of the strange distribution within a global 
fit of all parton densities. The wealth of neutrino and anti-neutrino data will 
also allow us to test the possible charge asymmetry of the strange sea [s 7^ s) 
predicted by some authors. 

In the present work the strong coupling constant is independently fixed to 
a value close to the world average. In a forthcoming paper we shall study the 
possibility of determining from the minimisation of the total of the fit, 
and discuss the stability of this determination and its correlation with the gluon 
density. 

The article is organised as follows. In section |^ we collect the main theoretical 
ingredients of the QCD analysis of DIS data. The re-evaluation of the neutrino 
and anti-neutrino differential cross sections is the content of section 1^. The fitting 



^In the GRV fit |Q the strange distribution is assumed to be zero at a very low scale 
and then entirely generated by the QCD evolution. 
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procedure is described in section ^. Finally, the results on cross sections, structure 
functions and parton distribution functions are presented in section |^. 



2 Deep inelastic scattering in QCD 



The differential cross section of neutral-current (NC) DIS of charged leptons {£), 
in the one-boson approximation and for moderate -C M|, is given by 



xy''F!-ix,Q'')+ 1-y 



2E J 



x,Q' 



dxdy 

If <^ M| charged-lepton NC DIS is essentially an electromagnetic reaction 
(that is dominated by one-photon exchange), the Z° contribution being totally 
negligible. In eq. ([^) = 1/137 is the electromagnetic coupling constant, E is 
the beam energy and Mjq is the nucleon mass. 

For charged-current (CC) neutrino (anti-neutrino) DIS one has 



dxdy 



G\M^E 
n 



xy'F,^^'^''{x,Q')+[l 



+ ( 



y 



\xE 



u(u)N . 



(2) 



where Mw is I^-boson mass and Gi? is the Fermi coupling constant. 

In QCD the structure functions F2, Fl = F2 — 2xFi and xF^ are given by the 
convolution of the parton distribution functions (pdf 's) with some perturbatively 
calculable coefficient functions. In the kinematic range covered by the analysis 
presented in this article the contribution of the b quark is negligible and the only 
heavy quark considered is charm. Therefore we split the structure functions in 
two components 



F,(x,g2)=F,,Kx,Q2)+F,,e(x,Q2 



2,3,L 



where Fj ;(x, Q^) is the light-parton contribution and Fi dXjQ"^) is the charm 
contribution. In the CC case the latter mixes charm with light quarks. 

The massive scheme that we adopt is the Fixed Flavour Scheme (FFS) |]1^ 
in which charm is a 'heavy' quark in absolute sense. This means that there is 
no such thing as the charm density function and charm is radiatively produced. 
Consequently the number of active flavours is set to 3 irrespective of Q^. The 
FFS has been shown to be more stable than the alternative massive scheme, the 
Variable Flavour Scheme (VFS) p3[, at moderate where most of the neutrino 
data lie Pl, ES. 
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F2,l 






x[|(u + u) + l{d + d + s + s)] 






x[^{d + d) + + M + s + 




up 


2x[u + d|Krf|2 + s|K,|2] 




un 


2x[d + u\Vud\^ + s\Vus\^] 




up 


2x[u + d\Vud? + slVusl"^] 


2mo\Vcd\^ + m\Vcs\^] 


un 


2x[d + nlKdP + 





Table 1: Leading-order expressions of -F2 {l- light sector; c: charm sector in the FFS). 
The X and arguments have been omitted. Only the slow-rescaled argument ^ has 
been explicitly indicated. 





xF^^l 


xF^^c 


up 


2x[d\Vud? + s\Vus? - u] 


ndmvcdi^ + 40\Vcs\^] 


un 


2x[u\Vud? + s\Vus? - d] 


2^[umvcd? + smvcsi^] 


up 


2x[u - d\Vud? - s\Vus?] 


-2mo\Vcd\^ + miVcsi^] 


un 


2x[d - u\Vud? - s\Vus?] 


-2mmcd\' + m\ycs\'] 



Table 2: Same as table ||, for xF^. 
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In the strong couphng, heavy quark thresholds are accounted for according 
to the prescription of Ref. ||2^ (with = 1.5 GeV and rrih = 4 GeV). 

The leading-order (LO) expressions for F2 and xF^ are collected in Tables |l] 
and 1^. Note that: 

• ^ = x{l + Q'^/ml) is the slow-rescaling variable; rric = 1.5 GeV is the charm 
mass. 

• Vij are the Cabibbo-Kobayashi-Maskawa matrix elements. We shall use 
\Vus\ = \Vcd\ = 0.224 and = = v^l - 

• The symbol (8> stands for convolution: 

dz 

f®9= —f{z)g{x/z) , 

J ax ^ 

where a = 1 for the light sector, a = 1 + for NC charm production, 

a = 1 + m'l/Q'^ for CC charm production. 

• The LO charm contribution to F2 in the neutral-current charged-lepton case 
is a 0{as) quantity in the Fixed Flavour Scheme. C^f^ is the LO Wilson 
coefficient for the photon-gluon fusion process Explicitly: 



Q') = I (g) C^2f\rnl/Q') ® xgif^') , (3) 

where the strong coupling is evaluated at the factorisation scale /i, and 
C2^g\z,m?^/ Q"^) can be found for instance in [[l7| . 

• The charm production is different in neutral and charged current DIS. In 
the former case it is given at LO by gluon splitting into a cc pair. In the 
latter case it is given by the direct process W^s c, with the slow-rescaling 
variable taking into account the effect of the charm mass. 

• At order the longitudinal structure function Fl = F2 — 2xFi vanishes. 

The QCD analysis performed in this paper is at the next-to- leading order 
(NLO) level (the renormalization scheme adopted is MS). NLO means 0{as) 
for the light sector and for the charm contribution to CC DIS, 0{al) for the 
charm contribution to NC DIS in the Fixed Flavour Scheme. Since Fl vanishes 
at order a^, for consistency with the treatment of the charm structure function, 
we include in our NLO analysis the order contribution to Fl, except for 
the strange-charm component of F^, for which the 0(af) longitudinal Wilson 
coefficients are not known yet. Anyway, this contribution has a very little effect 
in the kinematic domain of our analysis. 
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The hght-parton components of the structure functions have the form (for 
illustration we write only F2 for the electromagnetic case): 



f=g,g 



(4) 



The MS Wilson coefficients C^]^ and C^^J with i 
pendix I of Ref. [Q. As mentioned above, for Fl^i we consider also the 0{ 



2, 3, L can be found in Ap- 



contributions to the coefficient functions calculated in Ref. pOfPi 

The parton densities xf and xg are obtained by solving the DGLAP equa- 
tions at NLO [0. For the light sector we choose ^/Q^ as the factorisation and 
renormalisation scale. 

The NLO expression of the charm contribution to the NC structure functions 



IS 



+ 



+ 



a. 



(27r)2 ^ 9 

^ ' f=q,q 



.fD'',f{ml|Q')®xf{^,' 



(5) 



The NLO coefficients C"^'*-^-* and D'^'^^^ for and have been computed in 
the Fixed Flavour Scheme in For our calculations we have used the tables 



presented in [32 



The NLO strange-charm component of the structure functions is given by 
(omitting the Cabibbo-suppressed term and again writing only F2 for simplicity) 



2,c 



a. 



2{(l + ^i^2,f^ 



27r 



(6) 



The Wilson coefficients H'^'f^ and Hf^^^ 



2, 3, L) have been computed in p3 



and can be found, with a convention update, also in |3^, The factorisation and 
renormalisation scale for the charm structure functions is chosen to be a/Q^ + ml. 



3 Re-evaluation of neutrino cross sections 
3.1 Bin centre and radiative corrections 

We start with a general description of the procedure of bin centre and electroweak 
radiative corrections. We shall then give the details of the application of this 
procedure to the various data sets. 

■^We thank W. van Neerven for having provided us the code which computes the order aj. 
Wilson coefficients of F^. 
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The neutrino DIS cross-sections are determined experimentally in bins of three 
kinematic variables, say {x,y,E). The total cross section cr*°^ corresponding to 
the bin [xi,Xi+i], [yj,yj+i] and [^^,^^+1] is given by 



%+i fEk+i £^ 



<k = -^——dEdxdy 
I I I- dE dxdy ^ 



ijk 



c 



(7) 



where C is the number of scattering centres, Nij^ is the number of events corrected 
for detector effects and background contamination observed in the bin {i,j,k), 
d(j)/dE is the neutrino beam energy flux. 

In order to relate ajjl to the differential cross section, we invoke the average 
theorem: there exists at least one point {xi,yj,Ek) inside the bin {i,j,k) such 
that the following relation holds 



d'^a{xi,yj,Ek) ^ 0^ 



tot 

(8) 



where S.^k = XT^' Jy^^' Je^^' % dE dx dy is the bin surface. 

The differential cross section defined in eq. (§) must be corrected for elec- 
troweak radiation effects and translated to the bin centres (x^, yj, -E^). This is 
done by constructing an "experimental" Born differential cross section defined as 

d^'^expixl,yj,E^k) _ d'^a{xi,yj,Ek) 



dxdy dxdy 

S,,kd^~a^+''{xi,y'^^,El)/dxdy 



X 



iTKrrEr^'i^dEdxdy 



d'a^{xly^,El)/dxdy 
d^d^+^{xly'^,El)/dxdy ^ ' 

where the tilde symbol designates the quantities which are theoretically com- 
puted, and the upper-scripts B + R and B mean that these calculations are 
performed including or not, respectively, the higher-order electroweak correc- 
tions. The first term on the r.h.s. of eq. (P) is given by eq. (§). The second term 
embodies the bin centre corrections and requires a smooth parametrisation of the 
data. The third term incorporates the radiative corrections. 
Combining the last two terms terms we get 

d(^exp{xt.ypEl) ^Nijk^ dd^{xi,y''j,El)/dxdy 

dxdy ~ C ry^+^ r^^+^ ^^^^^l^piMdEdxdy' 

Jxi Jyj J E], dE dxdy «^ 



The electroweak corrections to charged current DIS were computed in [^, ^ 



In our analysis we used the program of Ref . , based on the results of |^ . At 
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the lowest order radiative corrections include the radiation of virtual and real 
photons from the charged lepton and quark legs, and the 7 — box diagram. 
We found that these corrections can reach ~ 20% in some kinematic domains. 
The correction factor - the third term in eq. @ - is a ratio, hence it is rather 
insensitive to QCD corrections . Thus, for simplicity, we computed it at leading 
order in QCD. 

We also evaluated the effects of the higher-order 7 radiation from the charged 



lepton leg [^, using the program HECTOR [0. This correction typically does 
not exceed ~ 0.5% and was applied only to the CDHSW data which are statisti- 
cally more significant. 

3.1.1 CDHSW 

The CDHSW Collaboration published |^ the vFe and vFe differential cross 
sections corrected for detector effects and background subtraction. The measure- 
ments are binned in x and y for 9 different values of the neutrino beam energies 
El between 23 GeV and 187.6 GeV. Neither the bin centre correction in x and y 
nor the radiative corrections were applied. Thus the full correction of eq. (|^) is 
needed. 

To evaluate the correction factors of (^) we cannot use any of the existing 
parton fits, as they do not account for nuclear effects. Thus we adopt an iterative 
procedure. In the first step, the bin centre correction and the radiative corrections 
are determined independently: the former by a parametrisation obtained in two 
different ways (see below); the latter by a standard fit (we use the LO GRV pdf's 



21]). Then, we perform a LO QCD fit to the corrected CDHSW cross sections. 



to the CDHS data and to the BCDMS, NMC and SLAC ill structure functionsQ. 



Using the pdf's of this fit we re-evaluate the full correction factor of eq. (y) and 
we iterate this step until the corrected differential cross sections get stable. In 
practice, to achieve the stability only two iterations are required. 

As we have mentioned, the smooth parametrisation of the data (required in 
the bin centre correction) is obtained by two different methods: by a fit of 
the CDHSW cross sections [fitting method), and ii) by an unfolding procedure 
[unfolding method). 

In the first method the published CDHSW uFe and uFe differential cross 
section data are fitted to 

.,+1 d''~a,(,.{x,y,E'; 



"'""^ ' ' dxdy, k = l,...9 



dxdy 



where 5"j^(p) is computed at LO using simple Buras-Gaemers-type pdf's |^ (which 
incorporate analytically the dependence). In parallel we adopted also the un- 
folding method (described in Appendix I). The difference of the results obtained 



Nuclear corrections to neutrino data are applied as explained in section 3.2 
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by the two procedures can be taken as a (partial) estimate of the uncertainty on 
the correction factor apphed to neutrino cross sections. We found that after two 
iterations the resuhs of the two methods are compatible within 1%. 

In fig. |l]a,b we show the total correction factors applied to the data of the 111 
GeV beam sample, as a function of yj for fixed x,. They vary between +6% and 
—4% and are roughly identical for neutrino and anti-neutrino beams. In fig. |l|c,d 
we show the contribution of the electroweak radiative corrections alone, which 
turns out to be the dominant one. 

As for the normalisation of the CDHSW data, in they were normalised 
using the average total cross sections of [Q, namely the ratios a"^*^ / Ey (a'^^'^/Ep) 
were assumed to be independent of the beam energy. This is a strong assumption 
in view of the fact that a linear rise with the energy is not experimentally excluded 
43]. It is therefore important to check the energy dependence^. 







4 


4 


(3 ± 5) • 10"^ 


0.98 ±0.05 


(-2.4 ±0.5) • IQ-^ 


1.00 ±0.04 



Table 3: Results of the minimisation of xt ^^'^ Xp- The hnear functions of the fits are 
defined as a/E =< a/E > {aiE + 02) where < a/E > is the average value (see text). 
Errors on the parameters are also given. 



To this end we performed a linear fit to the measurements (j^^^{Ejj) / E^ and 
a^^^{Ep)/ Ep of The results, renormalised by the average values 



< a^^yEy >= 0.703xlO-38cmVGeVand< (T^^VEp >= 0.331 xlQ-^^^cmV GeV, 
are given in Table 3 and shown in fig. ^ together with the one-standard-deviation 
band, computed according to the formulae of Appendix II (we call xl and xl 
the x^'s of these fits). A linear rise of a^^^{Ey) / Ey with E^^ is clearly compatible 
with the data. 

To take into account the effect of the uncertainty of this fit on the CDHSW 
data, we allow the parameters and to vary during the global pQCD fit by 
adding the term 

2 

« - ar)M^(aJ - aj) ± (af - af )M^(af - ap (11) 

«,i=i 

to the global x^ expression. Here a\ and a\ are the parameter values obtained 
from the preliminary linear fits and the matrices M,^- and M^j are defined as: 
M^. = {1/2) d\l/da\da^^, = (1/2) d\l/da^da^ . 

Another work done on the CDHSW data is the separation of correlated and 
uncorrelated systematic errors. The correlated systematic errors are dominated 

^We thank M.W. Krasny for having suggested to take the energy dependence of the CDHSW 
data into account in our analysis. 
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1^ by a possible shift of the hadronic energy by ±0.5 GeV. As no information 
has been pubhshed, we had to estimate the effects of this uncertainty on the cross- 
section measurements. In CDHSW, the beam energy E^y (Ep) is experimentally 
reconstructed via Ey = E^ + Ex — M^, where E^ is the measured outgoing 
muon energy and Ex is the measured hadronic energy in the final state. Thus 
a shift of the hadronic energy induces a variation of the kinematic variables: 
{x,y,Q^) — > {x±,y±,Q^). Our estimate of the relative differential cross section 
error 5^^^ induced by such shifts is therefore 



-(- Jxi Jyj dxdy 

ijk = p^+i p.+i ^ ^+J,y^E) dxdy 

Jxi Jyj dxdy 



where the cross sections a^^^ are computed as described above. The uncorrelated 
systematic errors are then obtained by subtracting quadratically the estimated 
correlated errors from the published systematic errors (J^f^^ 



The treatment of both types of errors in the fit is described in section ^. 

Finally, as CDHSW belongs to an old generation of experiments, one may 
worry about possible sources of errors at high y neglected in their analysis. In 
particular, one may question two points of the CDHSW analysis (see for more 
details): i) a. rejection cut against the dimuon events was used but no correction 
for it was applied; ii) the background coming from muon production in hadronic 
showers was not taken into account. Using the recent CCFR results on CTg^^/aJ'^'' 
W\ and on the production rate of muons in hadronic showers |H6[, we estimated 



the effect of these two sources of errors on the measurements and found it to be 
at the level of 1% at most. 



3.1.2 CDHS 

The CDHS Collaboration measured the v{i')H and z/(z/)Fe differential cross 
sections. Both the bin centre corrections and the radiative corrections were ap- 
plied. The latter, however, are incomplete. Only the 7 radiation from the muon 
leg was in fact considered. Hence our first step was to uncorrect back the pub- 
lished CDHS cross sections as follows 

d'^auncorr{xl,ylE^^) _ d'^ a pubiixj , y'j , E^ d^a^+^'{x^i,y^, E^)/dxdy 
dxdy dxdy d'^a^{x1,yp E'^ / dxdy 

where the upper-script B + fi means that only the radiation from the muon leg 
was included (to compute it we used a modified leading-order version of the 
HECTOR program). 
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Then we applied to (Pauncorr/dx dy the full radiative correction factor (i.e. 
the third term in eq. (§) 

d'^(Texp{x'i^ yh ^k) _ d^CTuncorrixt, d^d^jxl E^)/dxdy 

dxdy dxdy d'^a^^^{x1, y^, E^) / dxdy 

In the original paper Iron and Hydrogen data were normalised using 
ct'^-^V^^ = 0.625 X 10-38cmVGeV and a^^V^p = 0.3 x lO-^^cmVCeV. Since 
that time, the total uFe and uFe cross sections have been measured more pre- 
cisely . We therefore applied a new overall correction factor corresponding to 
a^^7E^ = 0.703 X 10-38cmVGeV and cr^-^V^p = 0.331 x 10-=^^cmVGeV [|3 



We have also improved the overall normalisation of the Hydrogen data combining 
the CDHS [| (j''(^)^/ct''(^)^'^ and BEBC ^^uim ^^uiu)Ne ^esuhs. Assuming no 



correlations between these experiments, we determined the absolute normalisa- 
tion of the CDHS Hydrogen data using: a"^ /E^ = 0.451 x lO^^^cm^/GeV and 
a"^ / Ey = 0.473 x 10~'^^cm^/GeV. The remaining error is 3.3% for the neutrino 
beam and 5.3% for the anti-neutrino beam. In summary, the uFe, uFe, vH and 
vH pubhshed CDHS data have been renormahsed by +12.5%, +10.3%, +14.5% 
and +20.6% respectively. 

3.1.3 BEBC measurements 

Among the BEBC publications only one ^ gives enough information to recon- 
struct the differential cross sections without any QCD assumptions. In this arti- 
cle, the corrected number of events is given in bins of x and y for a given range 
of Q^. As in the case of the CDHS data, the radiative corrections were applied 
incompletely, considering only the radiation from the muon leg. Hence we uncor- 
rected the published BEBC data analogously to what we did for CDHS (except 
that now we have number of events instead o differential cross sections) 



■\Tuncorr j^jpubl •^^i •^Vj -Emin dE dxdy ^ A\ 

- ^,k x.^, y^^, E ^ ''^y^y'"^^ dxd^ ' ^ ' 

Jxi Jyj JEmin dE dxdy 



Then we inserted N^Jl^"^"^ into eq. (|1^) and constructed the 'experimental' Born 
cross section 

d^^lM^y'-yEl) _ ^" d~a^{x'i,y^^,El)/dxdy 

dxdy C • p+i p.+i r^-^ d<t daB+nix,y,E) ^^^^^ ■ 

Jxi Jyj JEmin dE dxdy 



To compute the correction factors, the beam energy flux function of was 



usecl0 and the kinematic cuts applied in the data analysis were taken into account: 

^We thank U. Katz for having provided us the FORTRAN code computing the beam energy 
fluxes used in [Bsl. 
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g2 > 2GeV2, W'^ > 4GeV^ 5 < E < 160 GeV and > 5 GeV (E^ being the 
scattered muon energy). 

One should notice that: i) The beam energy flux is not known absolutely, 
hence the absolute normalisation of the data has to be determined by the global 
fit. ii) The systematic errors are not available. 

3.2 Nuclear effects 

Most of the neutrino DIS data come from experiments with nuclear targets 
(deuteron or heavy nuclei). Thus they have to be corrected for nuclear effects, 



which are known to be quite relevant (see for instance [4£, 



3.2.1 Deuteron 

We correct the calculated nucleon structure functions for nuclear binding, Fermi 
motion and off-shell effects using the results of the covariant approach of Mel- 
nitchouk et al. The neutrino structure functions F^^, and F^^ are 

treated analogously. 

The calculation of |51| describes only the high-x (x > 0.2) behaviour of the 



deuteron structure functions. At small x other mechanisms are at work (anti- 
shadowing and shadowing) but they are negligible in the x and regions of the 
deuteron data entering our analysis |]52| . 

3.2.2 Heavy nuclei 

All CDHSW and a large fraction of CDHS data are obtained from scattering off 
Iron nuclei. Since the theoretical understanding of nuclear effects in heavy nuclei 
is still uncertain and model dependent ^ , we adopt an empirical procedure 



to perform the nuclear corrections. 

The experimental i'{u)Fe differential cross sections are fitted to 

where -R^Jq ■* is the correction factor for the non-isoscalarity of Iron 

= 777T^ ,t ^,/7-, (A = 55.8, Z = 26) , (17) 

and da^^^^^^ is the 'isoscalar' Iron cross section which incorporates the nuclear 
corrections. Thus 



''We thank W. Mclnitchouk for having provided us the computer code of his calculation. 
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where Rn°ci is a function embodying the nuclear effects on an isoscalar target. 
^nud consists of two factors 

Rnud = Rpe/D ■ Riso ■ (18) 

The first factor is the Fe/ D structure function ratio 

Rpe/D = ^7d" ' (19) 

which is obtained from a fit to the pubhshed experimental data on F^'^/Fg^, un- 
corrected for isoscalarity. The second factor contains the isoscalarity corrections 



(Ff + Ff)/2 
[ZF^" + {A- Z)Fi^)/A 



Riso - ,^ . , . ^7T;7r77T • (20) 



Rpe/D is a function of x only because there is no experimental evidence of a 
significant dependence (for a recent study, see |5^). Theoretically, a higher- 
twist {i.e. power-like) dependence is expected at small x and ~ 1 GeV^ 
3^, but at larger Q^, in the region covered by our analysis, shadowing is a 



scaling phenomenon |56|, |5l 

The small-x {x < 0.1) uFe and i^Fe are excluded in our analysis. The reason 
is that in this region there are a number of uncertainty sources affecting the 
determination of Rnud- First of all, there is an unsolved discrepancy between 
the two measurements Rpe/D at small-x, namely between E665 ||5^ and NMC 



[Q. Second, the use of the charged-lepton DIS data to determine the Fe/D 



ratio of neutrino cross sections in eq. ([16|) is justified by the BEBC A^e target 



results being in good agreement with the NMC results on C (see |Jl|]), but 
the situation is experimentally not so clear for x < 0.1, where different nuclear 
corrections for charged lepton and neutrino structure functions are expected from 



a theoretical point of view [Q. Finally, the F2 and xF^ corrections might be 



different at small x Thus the cut at a; = 0.1 removes from our fits the more 
controversial kinematic region as for nuclear corrections. 

We have then performed a fit to the SLAC [|| BCDMS H Ff^/F^^ data 



in the range 0.1 < x < 0.65. For the function Rpe/Dix) we chose the empirical 
form 

Rpe/Dix) = «! + a2X + a^x'^ (21) 

where a^, i = 1, 2, 3 are the free parameters of the fit. The parameters min- 
imising the of this fit (we shall denote it Xfc/d) given in Table 4. The 
result of the fit is shown in fig. ^ together with the one-standard-deviation error 
band. 

In fig. ^a,b we plotted the isoscalarity and nuclear+isoscalarity corrections 
computed from the structure functions obtained with our fit (see section ^ . 
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1.091 ± 0.021 



Q2 

-0.34 ± 0.12 



as 

-0.09 ± 0.16 



Table 4: Results of the minimisation of Xpe/D (s^e text). 



4 The analysis 

4.1 Data entering the fit 

Our analysis includes the neutrino and anti-neutrino cross section data of BEBC 
[H, CDHS and CDHSW [0]. The number of collected events, the number 
of points and the kinematic domain covered by these experiments are listed in 
Table |. 

Besides the neutrino data, discussed at length in section |, the bulk of data 
entering our fits consist of structure functions from various charged-lepton DIS 
experiments: BCDMS HI pl, NMC Q (see Table 1). 

As for BCDMS, we do not use the merged structure function data, obtained 
by putting together measurements at different beam energies, but rather the data 
on the reduced cross section 



1 d^a 



A'koI^MnE Y dxdy 



MnXV 
E 



(22) 



for each beam energy (we retrieved the original data from ref. [34|). 





events 


Ax 


AQ2 


(GeV2) 


Ay 


CDHSW(z/Fe) 


640 000 


0.015 - 0.65 


0.2 


- 240 


0.037 - 0.87 


CDHSW(z?Fe) 


550 000 






n 


n 




2100 


n 




V 


V 


CDHS(z^i/) 


1100 


n 




V 


n 


BEBC{vD) 


12100 


0.03 - 0.65 


3 


- 75 


0.05-0.95 


BEQC{vD) 


5 400 










BEBC{vH) 


9 800 


0.03 - 0.65 


3 


- 75 


0.05-0.95 


BEBC{vH) 


4 900 


57 




)7 


17 



Table 5: Number of events and limits of the kinematic domain covered by the neutrino 
experiments considered in our analysis. 



To avoid a redundancy of data sets and, most of all, to limit the number of 
experimental parameters in the minimisation, the CDHS Iron data, which are 
in agreement with the CDHSW data, do not enter the fit. For the same reasons. 



the ZEUS F2 data are not included (ZEUS and HI 1994 data are consistent 
within 1 — 2%, as shown by a recent pQCD analysis 
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Ax 


AQ2 (GeV^) 


NMC (hp, fiD) 


0.008 - 0.5 


0.8 - 60 


BCDMS {^ip) 


0.07-0.75 


7.5 - 230 


Hl-94 [ep) 


3.2 • 10"^ - 0.32 


1.5 - 5000 


E605 (DY) 


0.12 - 0.4 


22.6 - 248 


NA51 (DY) 


0.18 


27.2 


E866 (DY) 


0.036 - 0.312 


30 - 164 



Table 6: Kinematic range of the charged-lepton DIS and DY data sets entering our 



Drell-Yan (DY) data are introduced in order to constrain the non-strange 
sea. Three measurements are used. We fitted the differential cross section for the 
reaction pCu —>■ fi'^fi'X measured by E605 [13|. The cross section is calculated 
at NLC)0. For Drell-Yan processes the charm contribution is small 
neglected in our analysis. No higher twist corrections are required 



and is 
39| since the 

kinematic domain covered by E605 avoids the phase space boundaries where these 
corrections are expected to be important|. 

The other Drell-Yan results that we use are the measurements of the DY 
asymmetry in pp and pD collisions by NA51 and by E866 [T^. These data 
constrain the ratio u/d. 

The kinematic cuts applied to the DIS data entering the fit are: 

• > 3.5 GeV2 and > 10 GeV^. 



In this region higher- twist effects are negligible |71|. Target mass corrections [72 



are also very small but have been taken into account in our calculations 
Due to the and cuts neither the E665 W^, nor the HI (1995) 



nor 



the SLAG |4T|, nor the ZEUS (1995) data l7| enter our fits. 

Finally, we reject the CDHSW data with x < 0.1. The reason for this cut 
is threefold: i) the systematic errors in the low-x region are large |Q; ii) the 



nuclear corrections at small x are not completely under control, as discussed in 
section |3l^ ; in) at low-x the CDHSW results disagree with the CCFR findings 
for the cross sections |76[ and for the structure functions 0. 



4.2 Fitting procedure 

The main steps of our fitting procedure are summarised below. For each iteration: 

^ We thank W. van Neerven for having provided us the code computing the order as Wilson 
coefRcients published in ||67|]. 



^This is not the case for the latest E772 Drell-Yan data 
consider them here. 



and that is why we do not 
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1. The pdf's are parametrised at a given value of the momentum transfer 
denoted Qq. We choose Ql = 4: GeV^. 



2. The DGLAP equations are solved numerically in the x-space (see |^ 
for a comparison of different NLO evolution codes). 

3. The evolved pdf's are convoluted with the Wilson coefficients to obtain the 
structure functions (see section 

4. Assuming that all experimental uncertainties are normally distributed the 

is computed 



2 \ ^ \ ^ L exp 



exp dat ^dat,stat ~^ ^ dat,uncor 



exp\2 



exp exp k 

2 

+ ^[(ar-ar)iW^K-aP + (<-af)iW^K-«f)]' (23) 

where O stands for the observables (structure functions and differential 
cross sections). The first two sums run over the data [dat) of the various 
experiments {exp); a exp is the relative overall normalisation uncertainty; 
o'dat,stat and <Jdat,uncor are the statistical error and the uncorrelated system- 
atic error, respectively, corresponding to the datum dat; Vexp is the number 
of standard deviations corresponding to the overall normalisation of the 
experimental sample exp; ^k^^i^s^k^) is the relative shift of the datum dat 
induced by a change by s^^^ standard deviations of the k*^ correlated sys- 
tematic uncertainty source of the experiment exp. It is estimated by 

C^tdatf exp _ ,1 N _ (^datf exp _ _i \ 

s^datt exp\ _ ^exp\^k ~^ ) ^expK^k ^ exp , 

y^k ) - 2ndat + 

^exp 



2Qdat 

^ exp 



^exp\2 



where O'^'^^^^s^^^ = ±1) is the experimental determination of Of"* obtained 
varying by ilo" the k^^ source of uncertainty. The last term in eq. (^3]) has 
already been discussed in section |3.1.1| (eq. (pA])). 

Notice that even though the parameters Uexp, s^k^, a^ and a^ are obtained 
from the global minimisation, they do not enter in the counting of the 
degrees of freedom since they are determined by the counter-terms. 

The correlated systematic uncertainties are taken into account whenever the 
information about them is available. This is the case for HI [Ol, BCDMS 
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33 and NMC [ITOl. For CDHS and BEBC no information is available and for 



E605 |T3[ the uncorrelated systematic uncertainties are dominating. The 
systematic uncertainties of this later experiment has then been added in 
quadrature and included in cr^^^ ^^^^^ of eq. (|23|) . As already pointed out in 
section |3.1.3| , the systematic uncertainties of the BEBC data are not known. 
Since these data cover the same region as the CDHS data, for both data 
sets we have taken into account only the statistical uncertainties (this has 
a negligible effect on the minimisation, since the statistical significance 
of these samples is rather small). 



5. The MIGRAD algorithm of the MINUIT program [|7^] is used to minimise 
the x^- 



Given the importance of nuclear effects in the treatment of the neutrino data, 
in parallel to the main fit, we performed, as a check, another fit, in which the 
nuclear parameters of eq. ( PT| ) are not constrained to the values di obtained by 
the independent parametrisation of the Fe/D structure function ratio described 
in section |3.2.2| , but are readjusted a posteriori. This is done by adding the 
counter-term 



E 



{oi - ai)Mf^'^{a.j - aj) 



(24) 



to the expression of the global pQCD analysis. In eq. (^ M^'^^ = (1/2) dW^j^/ daidaj. 
We found that the two fits give very similar results. 



4.3 The parametrisation 

Imposing the isospin symmetry leads to the following relations among the pdf 's: 



P = d'' = M, dP 



u 



rf, = rf" = M, dP 



U, S'' 



In our main fit, that we call f itl, the pdf's 



u 



■■ S = S, S'' 



= s = s. 

d, M, d, s, s 



and g (the gluon density) are parametrised at Qo = 4 GeV^ as follows: 



xu^{x, = Au^ x^"" (1 - x)^"- (1 + Du, x^"") 
xd,{x, Ql) = Aa, x""^^ (1 - xf''^ (1 + Dd, x^'^-) 
x(m + rf)(x, Ql) = A+ x^+ (1 - xf+ (1 + D+ x^+) 
x{d-u){x,Ql) = A_ x^- {l-xf' (1 + L>„ x) 
xs(x, Ql) = xs{x, Ql) = As x^' (1 - x)^= (1 + D,, x^") 



xg{x,Q. 



AgX^^ (1 - X 



« 'l + Dgx""'] 



(25) 
(26) 
(27) 
(28) 
(29) 
(30) 



This parametrisation form is similar to that used in (we refer to this article 
for a justification of this choice). 
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Generally it is also assumed s = s. The data samples used in the existing 
global analyses [|l^, ^ cannot resolve in fact s and s independently. In our 
case, the information coming from neutrino and anti-neutrino differential cross 
sections allows testing the hypothesis s = s. We thus performed another fit, 
called f it2, allowing for a charge asymmetry in the strange sea, s s. 

Some of the parameters in eqs. (p5|-^) are determined by physical constraints. 
One normalisation factor, say Ag, is fixed by the momentum sum rule, ^^{xg + 
+ qi))dx = 1. The two normalisation parameters Ay_^ and A^^ are fixed 
by the number sum rules Uydx = 2 and d^dx = 1. 

While the intermediate-x and large-x shape of the strange distribution is well 
constrained by the data entering the fit, the small-x behaviour is not. Thus we 
set Bs = We also set B^^ = Bd^ = as suggested by Regge theory. It is 
to be mentioned that, from a statistical point of view, these constraints do not 
worsen the cind are also required in order to get an invertible second derivative 

matrix. 



5 Results 

We present now the results of our fits. The strong coupling is fixed at the value 
as{M^) = 0.120 which is close to the world average pT] . 



The contributions of the different data sample to the s-^e given in table ^ 
The total per degree of freedom is excellent for both fits. The amount of 
systematic corrections for each datum, that is the value of the term [I'expO'exp + 
Eik *(sr^)] in eq. (|2|), is between -6% and +12% for all fits. 

The parameters of eqs. (]25|-pOD are listed in table || and the parton densities 
are shown in fig. ^, where they are compared to the results of the other global 
fits. The momentum fractions of the various partons are listed in table |^. 

The parton distributions in fig. |^ are accompanied by the error bands com- 
puted as explained in Appendix II. These do not take into account the uncer- 
tainties related to the functional choice of the pdf's, nor those inherent to the 
treatment of the errors, which are assumed to be normally distributed. The 
meaning of the error bands of our pdf's is the following. Once a specific form for 



the pdf's is chosen and the constraints described in section |4.3| are imposed, the 
error bands correspond to an increase of the by one unit. Thus their width 
is determined not only by the abundance and the precision of the data but also 
by the constraints on the pdf's. This explains why the error bands may be small 
even in kinematic regions where there are no data. 

With respect to the other parametrisations, our gluon density turns out to 
be higher at intermediate x. This can be due to the fact that we do not use the 
prompt photon data, which tend to favour a larger glue at high x and a smaller 
one at intermediate x. These data are still quite controversial and some of them 
seem to be in disagreement with the QCD predictions (for a recent discussion see 
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Experiments 


#pts 






# exp 






fitl 


fit2 


param 


BCDMS(IOO) cr^P 


94 


108.0 


110.8 


12 


BCDMS(120) cr^P 


99 


81.6 


81.0 




BCDMS(200) a^'P 


79 


92.7 


90.7 




BCDMS(280) at'P 


76 


89.7 


87.6 




BCDMS(120) a^'^ 


99 


96.9 


93.1 




BCDMS(200) 0-^-^ 


79 


93.5 


88.8 




BCDMS(280) 0-^^ 


76 


64.6 


63.7 




BEBC a'^P 


68 


65.2 


67.6 





BEBC a^P 


49 


76.4 


76.5 




BEBC a"^ 


70 


65.3 


65.4 




BEBC a^^ 


49 


49.7 


46.7 




CDHS a'^*' 


45 


50.6 


48.6 


2 


CDHS a^P 


42 


53.7 


53.2 




CDHSW fj^^^ 


494 


264.9 


250.4 


5 


CDHSW cr^-^^ 


492 


274.2 


277.5 




E605(DY) 


136 


104.9 


104.1 


1 


E866 ^Ipy 


11 


8.2 


8.2 


1 


rli(^y4-SVXJ 




OP. 7 


ZD. ( 


7 


B.l(9i-nvx) F^^ 


156 


180.6 


180.3 




NMC(90) ^2''^ 
NMC(120) Fs''^ 


34 


32.6 


32.5 


16 


46 


67.0 


67.3 




NMC(200) Fs''^ 


61 


99.4 


99.3 




NMC(280) F2^^ 


68 


106.1 


105.6 




NMC(90) F^^ 


34 


24.0 


23.8 




NMC(120) F^^ 


46 


50.6 


50.8 




NMC(200) F^'^ 


61 


62.6 


62.7 




NMC(280) F^^ 


68 


94.8 


95.6 




NA51 Ady 


1 


2.7 


2.7 





Total 


2657 


2430.8 


2405.0 


44 



Table 7: Contribution to the global , and number of points, of the data samples 
entering the fits. The values of the individual x^'s do not include the normalisation 
and correlated systematic shifts. The last column indicates the number of experimental 
parameters: overall normalisation of the data sets and correlated systematics. The 
contributions of the experimental parameters to the amount to 43.6 for fitl and 
40.5 for f it2. 
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fitl 


fit2 


y 


7.72 


7.3 


y 


0.089 


0.081 


Co 

y 


21.26 


20.11 


y 


12072 


6990 


y 


4.17 


4.0 


LLv 


1.43 


2.05 


LLy 


0.49 


0.55 


llv 


3.60 


3.75 




4.47 


3.59 


En, 


0.81 


0.97 


Ad 

Liu 


1.02 


1.40 


Bd 


0.49 


0.55 


Cd 

Liu 


6.03 


6.63 


Dd, 


23.06 


32.63 


Ed. 


1.76 


2.07 


A+ 


0.071 


0.075 


B+ 


-0.245 


-0.240 


c+ 


8.31 


8.62 


D+ 


11.29 


13.30 


E+ 


0.88 


0.97 


A- 


0.11 


0.12 


B_ 


0.49 


0.55 


C- 


16.08 


16.31 


D_ 


-55.62 


-58.46 




1 


1 


As 


0.064 


0.066 


Bs 


-0.245 


-0.240 


Cs 


5.31 


5.59 


Ds 


443 


11354 


Es 


8.26 


12.04 


As 


0.64 


0.066 


B-s 


-0.245 


-0.240 


Cs 


5.31 


5.44 


Ds 


443 


339 


Es 


8.26 


7.39 



Table 8: Values of the parameters of the pdf's at Qo = ^ GeV^. 
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9 




dv 


u 


d 


s 


s 


5 

GeV2 


42.9±0.4 
42.8±0.5 


27.9±0.3 
27.9±0.3 


11.2±0.2 
11.2±0.2 


3.0±0.1 
2.9±0.1 


3.9±0.1 
3.8±0.2 


2.1±0.2 
2.3±0.2 


2.1±0.2 
2.2±0.2 


20 
GeV2 


46.6±0.3 
46.5±0.4 


24.6±0.3 
24.7±0.2 


9.8±0.2 
9.9±0.2 


3.2±0.1 
3.1±0.1 


4.0±0.1 
3.9±0.1 


2.4±0.2 
2.6±0.2 


2.4±0.2 
2.4±0.2 


100 

GeV2 


49.4±0.3 
49.3±0.3 


22.1±0.3 
22.1±0.2 


8.8±0.2 
8.8±0.2 


3.3±0.1 
3.3±0.1 


4.0±0.1 
4.0±0.1 


2.6±0.2 
2.8±0.2 


2.6±0.2 
2.7±0.2 



Table 9: Fraction of the total nucleon momentum carried by the partons for three 
values of Q^. The results of f itl (upper row) and f it2 (lower row) are displayed. The 
errors are computed as explained in Appendix II. 



B^). Also their compatibility with the DIS data is still an unsettled issue. Our 
glue is determined by DIS measurements only and might be in slight disagreement 
with the prompt photon data at medium and high x. We will explore this problem 
in a future work where we will make use of an enlarged dataset. 

We also find a shght discrepancy in the distribution between our results 
and the other global fits. This is not surprising, since in our parametrisation we 
did not include the CDF data on the lepton asymmetry in W production at the 



Tevatron |83|. These data constrain the u/d ratio at a: = 0.05 — 0.1, which is 
precisely the region where some difference can be seen between our dy{x) and the 
CTEQ and MRST distributions. 

The structure function F2 measured in different charged lepton DIS exper- 
iments is shown in figs. ^ |^ and || with the curves of f itl. Our fits for the 
Drell-Yan data are presented in figs. ^ and |10|. An excellent overall agreement is 
observed. In fig. |TO|.b we plotted the ratio d/u at = 30GeV^, together with 
its error band. Notice that the fit, which is dominated by the E866 data, yields 
a. d/u ratio which lies below the NA51 determination. 

The reduced v (z/) differential cross section 

c?V 27r(M^ + g^)^ 1 2 M^xy 

d^y= GIM.M^^E Yd^y^^''''^ = ^^^^-y^ " ' ^'^^ 

computed from f itl is compared to the data in fig. [Tl| (Hydrogen target), fig. |12 
(Deuterium) and fig. |13| (Iron). The BEBC and CDHS Hydrogen and Deuterium 
data are well fitted down to small x. These data contribute to constrain the 
valence distributions, without being affected by nuclear effects. Fig. ^ also 
shows that the CDHS Iron data, though not entering the fit, are well described 
by it. As for the CDHSW rejected data (x < 0.1), discrepancies with the fit 
appear only in the first bin x = 0.045 of figs. |l^. 

Figs. |TB|, H, ^ and |^ and the results of table |^ show explicitly that the u and u 
Iron data are compatible with the data on F2 coming from NC charged-lepton DIS 
even in the region 0.045 < x < 0.2. By contrast, in ^ a sizeable discrepancy 
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is found between the CCFR F2 and the NMC F2 . From the analysis performed 
here, which is done - we recall - on neutrino cross sections, no disagreement 
emerges between charged-lepton and neutrino DIS measurements. We checked 
that including in our analysis the CCFR structure functions the fit worsens. 
We also found that the F3 data are much better described than the F2 data. 
We do not give any quantitative information about the of this particular fit 
since we considered the statistical errors only (the use of the total systematic 
uncertainties is not recommended by the CCFR Collaboration). In conclusion, 
some incompatibility seems to exist between the CCFR structure functions and 
all other charged-lepton and neutrino DIS data. A conclusive word on this matter 
could come from the analysis of the CCFR cross sections, which are unavailable 
at the moment. 

The longitudinal structure function Fi and the longitudinal to transverse 
cross section ratio R 



r2 ^2 



computed using the pdf's of fitl (we remind that all our fits include target 
mass corrections) are compared to the NMC |]T0| and BCDMS |1T[ results in 



fig. 0a, to the HI measurement in fig. |l^.b and to the CDHSW measurement 



[0 in fig. |T^.c. Again, one can see a good agreement between our fit and the 
experimental results. 

The beam energy dependence of the CDHSW uFe and uFe data, resulting 
from the global minimisation, is shown in fig. ^ The deviation from the 
independent linear fit of section is small for uFe data and larger for uFe. 

5.1 The strange sea density 

Let us concentrate now on the strange sea density. Due to the lack of data able 
to constrain it, this distribution plays a lesser role in the existing global fits. 
Two of them (CTEQ ^ and MRST |2|), guided by the results of the CCFR 



determination of s{x) ||2^, in particular by the CCFR value of the strange-to- 
non-strange momentum ratio k = {x{s + s)) / {x{u + d)) ^ 0.5, impose 

s{x) + s{x) = ^ [u{x) + d{x)] . (32) 



In the GRV analysis [gl|], instead, the strange distribution is set to zero at 
the input scale, and then radiatively generated. 

The abundance of our neutrino and anti-neutrino data sets allows us to fit 
s{x) with no extra constraints. The resulting s distribution of fitl is shown, 
with its error band, in fig. ^, where it is compared with the other fits. Notice 
that it turns out to be closer to MRST and CTEQ at large x, and to GRV at 
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low X, where it is however much less constrained. From table ^ one sees that in 
f itl K = 0.67 at g2 = 20 GeV^. 

In fig. |15| s{x) is plotted at three different values, together with the CCFR 



results p2[. The agreement is good, although we must recall that the CCFR 



extraction of the strange distribution from dimuon production in iy(iy)Fe DIS 



has been criticised in many respects JTS], ^ |8^ 



Fig. |16| shows quantitatively how f itl is favoured with respect to another 
fit, that we call f itlb, in which the constraint (32) is imposed. The total of 



fit lb is 2492.4, higher than the of f itl by 62 units. In fig. [T^ we plotted the 
mean value of the so-called pull distribution, as a function of x, for the CDHSW 
data: 

<P'^^^>-da.= ^^-r—^2^ , ^ = > (33) 



dat,stat ' dat,uncor 

where Nexp{xdat) is the number of data at x = Xdat for the experiment exp. In our 



case exp stands for CDHSW(z/Fe) or CDHSW(z/Fe). One can see from figs. \16 



a 



and |T6|.b that it is the u data which tend to favour f itl with respect to f itlb. 

A very interesting question is whether the strange distribution is equal or not 
to the anti-strange one. The usual assumption s{x) = s{x) is not dictated in fact 
by first principles. 

We thus looked for a possible charge asymmetry of the strange sea performing 
a fit, called f it2, in which we release the constraint s = s. From table |^, one 
can see that the attempt of disentangling the s and s distributions is justified by 
the abundance of anti-neutrino events in our data set. 

The parametrisation (^) for s = s is replaced in fit 2 by two independent 
functions for s and s 

xs{x, Ql) = A, (1 - xf^ (1 + x""'), (34) 
xs{x, Ql) = AsX^' (1 - xf' (1 + Dsx'^'). (35) 

We set As = Ag and Bg = Bg, and we fix one more parameter by imposing 
{s — s) dx = (no net strangeness). 
The main results of f it2 are: 

• The minimum decreases by 25 units with respect to f itl (see table 0). 
Hence the choice s 7^ s is slightly favoured. 

• The strange distribution turns out to be harder than the anti-strange one. 



The difference s — s is shown in fig. 17. a with its error band. In fig. 17. b 



we plot the ratio s/s at Q"^ = 20 GeV^. 

The momentum fractions of the fit2 partons at different values are 
listed in table |^. The momentum fraction (xs) is larger than (xs). No- 
tice also that f it2 favours a higher value for the strange-to-non-strange 
momentum ratio k, with respect to the fit with s = s. 
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Let us comment now on a previous test of the strange sea asymmetry. In 
2^ the CCFR Collaboration found no evidence for s ^ s. In their analysis the 
constraint s{x) = [u{x)+d{x)]xA{l—x)'" was imposed, which limits the flexibility 
of the fit. Moreover, the dimuon sample of CCFR does not cover the high-x region 
and consists only of 5000 neutrino events and 1000 anti-neutrino events. Our 
analysis does not have these limitations (our data sample is much more balanced 
between neutrino and anti-neutrino events and no extra constraints are set but 
those ensuring no net strangeness), and gives a more precise result on s/s, as one 



can see from fig. |17 



In order to understand how the difference s — s is constrained by the data, 
let us consider the quantity 



AnxjM^ + Q 



2\2 



(36) 



dxdQ"^ dxdQ'^ 

The flavour content of A'^"'^ is more evident in the parton model where it reads 

A'^'" oc xs{x) — xs{x) + Y_[xu^{x) + xdy{x)] , (37) 

with y_ = 1 — (1 — 1/)^. The v — v cross section difference ( ^6]) is plotted as a 
function of YL, at fixed x and Q^, in fig. |18|. Comparing f itl and f it2, one can 
see that their results deviate with increasing x (being very close to each other for 
X < 0.3). At high-x f itl undershoots the CDHSW values of A^ for all bins. 
Fig. ITBI shows that the CDHSW data are more precise at high-y. Hence it is this 



region which drives the result on s — s. Looking at fig. [T^.c,d one sees also that 
it is the V data which favour f it2 with respect to f itl. 

Since the nuclear corrections applied to the CDHSW data are sizable (see 
fig. §.b), one may naturally ask to what extent our results are affected by the 
uncertainties on the evaluation of these corrections. The Fe/D ratio, as we 
have seen (see fig. |] and the check described at the end of section ^]2|), is well 
determined by the charged lepton data. Moreover, it factorises in the v — u 
difference. Hence this component of the nuclear correction is rather harmless. 
The isoscalarity corrections are instead different for neutrino and anti-neutrino 
observables and quite large. Looking at eqs. (|T7| , |20|) one sees that the isoscalarity 
ratio cannot exceed A/{2Z) = 1.073 for Iron. Fig. ^a shows that the maximal 
value we obtained is not far from this upper bound. We checked that reducing the 
isoscalarity ratio the resulting s — s difference gets larger. Thus we conclude that 
our results on the strange sea asymmetry are not spoiled, at least qualitatively, 
by the uncertainties on the isoscalarity corrections. 

Theoretically, a charge asymmetric sea is accounted for by introducing a dis- 



tinction between extrinsic and intrinsic qq pairs |87]. The extrinsic sea con- 
sists of short-lived quarks and anti-quarks produced by QCD hard subprocesses 
(bremsstrahlung and gluon splitting). It is evident that the extrinsic component 
of the sea cannot be charge asymmetric. On the other hand, the intrinsic qq pairs 
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exist over a longer time scale and are associated with nonperturbative phenom- 
ena. These pairs are still produced by gluon fragmentation but have time, before 
recombining, to interact with other partons. They represent higher Fock states of 
the nucleon {\qqqqq ■ ■ ■)) and manifest themselves in meson-baryon fluctuations. 

There is no fundamental principle forbidding a possible charge asymmetry 
of the intrinsic sea. Actually, there are reasons to believe that such asymmetry 
should indeed be a property of the strange and charmed sea [^, If the 

strange (or charmed) sea is asymmetric at low due to some nonperturbative 
mechanism, the QCD evolution simply preserves this asymmetry because s — s (or 
c — c) evolves like a non-singlet distribution and its first moment is constant. An 
interesting feature of the intrinsic sea is that it tends to exist at relatively large 
values of x [^, |8^, corresponding to the most energetically favoured configuration 
of the nucleon light-cone wave function. 

In the simplest model 



attributed to the p AK 
mesons have relatively small masses 



the production of the intrinsic strange sea is 
fluctuation. Due to chiral symmetry pseudoscalar 

|58| , the average x of 



As a consequence 

the s anti-quark in the K is smaller that the average x of the s quark coming 
from the A. Thus the s distribution is expected, on quite general grounds, to 
be harder than the s distribution. This expectation has been substantiated by 
explicit calculations in ||89|, Q (for other models see ||9l|]). Our results on the 
strange and anti-strange sea density are, at least qualitatively, in agreement with 
the predictions of the intrinsic sea theory. 



6 Conclusions 

We have presented a global next-to-leading order QCD analysis of a large set 
of DIS data, including the (properly re-evaluated) neutrino and anti-neutrino 
differential cross sections of BEBC, CDHS, CDHSW, the charged-lepton structure 
functions of NMC, BCDMS and HI, and the Drell-Yan data of E605, NA51 and 
E866. The full use of the information on the nucleon structure embodied in 
neutrino DIS observables and a proper treatment of the experimental systematic 
uncertainties are the main novelties of our approach. In particular, the large- 
statistics CDHSW Iron data allow disentangling the strange sector from the non 
strange one: this leads to a consistent determination of s{x) within a global fit, 
similarly to what happens for the other parton distributions. 

The charm mass effects, whose relevance for an accurate determination of the 
parton densities is a recent firm acquisition, have been consistently treated in a 
massive factorisation scheme, the Fixed Flavour Scheme. 

We found no evidence of any discrepancy between the neutrino data we con- 
sidered in our fit and the charged-lepton data. A complete and unambiguous 
analysis of all neutrino experiments and a conclusive check of the compatibility 
of their data with the charged-lepton data would be possible only if we could use 
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the CCFR differential cross sections, which are unfortunately not available. 

The large statistics of anti-neutrino events in our data sample allowed us to 
test the hypothesis of a charge asymmetric strange sea: s ^ s. We found some 
evidence for such an asymmetry. The qualitative features of the resulting s and 
s distributions (namely, a large x tail at low and s{x) softer than s{x)) agree 
with the expectations of the intrinsic sea theory. 

Finally, we outline some developments. As for the strong coupling, in the 
present work we took a pragmatic attitude, using the world average value. A more 
systematic study, in which we will extract from the DIS data and investigate 
its correlation with the gluon density, will be presented in a forthcoming paper. 

Another important development that we have in mind is the comparison of 
the results obtained in the two QCD massive schemes, FFS and VFS. 

Finally, one should envisage some independent, and perhaps more direct, ex- 
perimental test of the charge asymmetry of the strange sea for which we presented 
here only an indirect statistical evidence. 
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Appendix I: Unfolding procedure 

A two dimension net in x and y is defined. The nodes correspond to the cen- 
tres of the experimental bins. Following [^, we define a basis of continuous 
functions 0j(a;) and 'ipjiy) such that 0i(a;^) = du and ipjiVm) = ^jm (where Su is 
the Kronecker symbol). Provided these functions are determined using a spline 
interpolation the differential cross section can be written as 

^ 1=1 m=l ^ 

Applying the average theorem, eq. ( pSf ) leads to 

1=1 m=l 

with 



'^ijk ^ ^ ^ ' y 
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where Nijk is the number of events experimentally observed and C is the number 
of scattering centres in the target. 

The differential cross sections at the bin centres (fa'^^'^^^^{x'[,y^,El)/dxdy 
are obtained by reversing this system of Hx x Uy equations. Numerically, the 
more stable results were obtained using the first order spline. In this case, (j)i{x) 



and ipj{y) are the Lagrange (or 'hat') functions ||92 



= ^ {xi+i - x)l (xi+1 -Xi) ,Xi<x< Xj+i (40) 
, otherwise. 

The same expression holds for ipj{y). 

However, using eq. (POf), one can see that the system of eq. (|39| ) is incomplete 
if Xi 7^ or Xn^ 7^ 1 or ?/i 7^ or ?/„^ 7^ 1. Indeed, the boundary bins receive some 
contributions (in the r.h.s. of eq. (PU])) coming from the bin neighbours which are 
not included in the measurements. To overcome this difficulty, we have defined 
artificial extra bins [a;„^,x„^+i] Vy, [yo,yi] Vx and [yny,yny+i] Vx. Then we fixed 
the cross section at the centre of these new bins to a certain fraction A of the 
closest cross section measurement. Under this modification, eq. (|39D becomes an 
inhomogeneous system 

'"■y Tlx \ 

m — 1 / — 1 / 

j2 u{u)Fe/ c c rpc\ 



m=l 1=1 

Tlx "TLy 

^ ^ y dxdy 

1=1 m=l ^ 



and the differential cross sections at the bin centres are obtained by inverting this 
system with A = 1. 

In order to test the sensitivity of our results to the choice of A, we have set A 
to some extreme values: A = 2, A = 1/2. We observed a variation of the order 
of ~ 20% in the highest x bin for all y bins. But we point out that the results 
are completely stable in the first x < 0.1 bins for all y bins. This method is 
then used to perform the definitive bin centre correction in these particular bins. 



as mentioned in section |3.1.1| . The uncertainty of the method is estimated by 
repeating the procedure using a second order spline interpolation to define the 
basis functions 0i(x) and ipjiy)- The ratio of the two bin centre correction factors 
is then taken, bin by bin, as an estimate of the uncertainty due to the method, 
and added in quadrature to the uncorrelated systematic error of the measurement 
(it is of the order of < 2% and reaches 7.5% at high y). 
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Appendix II: Error bands 

We describe here the formula used to calculate the error bands shown throughout 
the paper. Calling p = {pi, ...,Pn} the vector of the free parameters of the fit, 
the error band of a given function / is given at each (x, Q"^) point by p3[ : 



Af{x, Q'; po) = \fix, Q'; p, + Ap(x, Q')) - f{x, Q'; p, - A,{x, Q'))\ (42) 

where p^ denotes the parameter set minimising the cind the vector 
Ap{x,Q^) = {Ap,(x,Q2)^...^Ap„(x,Q2)} is given by 

A (X Q^) = M-^OJ{x,Q^;p) 

' ^d,f{x,Q^-p)M-^d,f{x,Q^-p) 

with Mij = (1/2) d'^x^/dpidpj and dp = {d/dpi, ...,d/dpn}- 
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Figure 1: a) Correction factor applied to CDHSW uFe data of the 111 GeV 
beam energy sample, as a function of y for various x bins; b) same as a) 
but for CDHSW uFe data; c) electrowcak radiative correction factor Srad = 
{dPa^^^l dxdy) / {cPa^ / dxdy) — 1 in percent for CDHSW uFe; d) same as c) but 
for CDHSW uFe. 
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Figure 2: Total cross sections of uFe (top) and uFe (bottom) from CDHSW. 
The shaded areas are the one-standard-deviation error bands corresponding to 
the hnear fits described in section |3.1.1| . The curves are the results of f itl. The 
error bars correspond to the quadratic sum of statistical and systematic errors. 
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Figure 3: Fe/D structure function ratio. The full line is the result of a sec- 
ond order polynomial fit and the shaded area is the corresponding one-standard- 
deviation error band (see section ^]2|). The error bars correspond to the quadratic 
sum of statistical and systematic errors. 
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Figure 4: Corrections applied to Iron target data: a) isoscalarity corrections; b) 
nuclear and isoscalarity corrections. 



37 




Figure 5: The parton distribution functions of fitl, with their error bands, 
compared to the GRV98 (sohd hue), CTEQ4M (dashed hue) and MRST (dotted 
hne) fits. The results of fitl for u^, u and J - hardly visible - nearly coincide 
with the MRST and CTEQ4M results. 
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Figure 6: HI data vs. f itl. BCDMS and NMC Hydrogen target data belonging 
to the X domain covered by HI are also shown. These two data sets have been 
rebinned into the HI x bins for plotting purpose. The data are renormalised 
by the overall normalisation factor determined by the fit and they arc plotted 
with an additive bin constant c{xi) = 0.6 * {i — 0.5) corresponding to Xi = 
{0.32, ...,8.10""^}. The error bars represent the quadratic sum of statistical and 
uncorrelated systematic errors. 
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Figure 7: NMC Deuterium target data vs. f itl. The different beam energy sam- 
ples are shown separately. The data are renormalised by the overall normalisation 
factor determined by the fit and multiphed by a constant c{xi) — {4.8, 4, 3.2, 2.5, 
2, 1.5, 1.2, 1, 7.5, 5.2, 3.7, 2.5, 1.7, 1.2, 1, 1} corresponding to Xi = {0.0125,..., 
0.07, 0.09,..., 0.5}. The data have been re-binned into these Xi bins for plotting 
purpose. The error bars represent the quadratic sum of statistical and uncorre- 
lated systematic errors. 
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Figure 8: BCDMS hydrogen target data compared to f itl. The different beam 
energy samples are shown separately. The data are renormalised by the overall 
normalisation factors and the correlated systematic shifts determined by the fit. 
The error bars correspond to the quadratic sum of statistical and uncorrelated 
systematic errors. As we use the BCDMS differential cross-section data in our 
fits, the F2 data shown in this plot have been determined by using Fi from f itl. 
The data are renormalised by the overall normalisation factor determined by the 
fit and multiplied by a constant c{xi) = {3, 2.5, 2, 1.5, 1.2, 1, 1, 1, 1, 1, 1} 
corresponding to Xi = {0.07,..., 0.275, 0.35,..., 0.75}. 
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Figure 9: E605 data vs. f itl. The data have been renormahsed by +8%, as de- 
termined by f itl. The data are muhiphed by a constant c{xFi) = {10^, 10^°} 
corresponding to Xpi = {0.275, —0.125}. The error bars represent the 
quadratic sum of statistical and systematic errors. 
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Figure 10: a)E866 vs. fitl. The data have been renormahsed by +2% as 
determined in fitl. The error bars represent the quadratic sum of statistical 
and systematic errors. Notices that the data points correspond to different values 
of Q^. The line are thus an interpolation between the calculations of fitl. b) 
NASl's result for u/d compared to fitl and various parametrisations (CTEQ4M 
does not include the E866 data). 
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Figure 11: CDHS and BEBC Hydrogen data vs. f itl. The BEBC data have 
been rebinned into the CDHS bins for plotting purpose. The data are renor- 
mahsed by the overall normalisation factor determined by the fit and they are 
plotted with an additive bin constant c{xi) = {11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1.5, 1, 0.5} 
corresponding to Xi — {0.03, 0.65}. The error bars represent the statistical er- 
rors. 
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Figure 12: Same as fig. [Tl| but for BEBC Deuterium data. Here the additive 
bin constants are c(xj) = {8.5, 7, 6, 5, 4, 3, 2, 1.5, 1, 0.5}. 
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Figure 13: Same as fig. |TT] but for tlie CDHSW 65.7 GeV beam data sample. 
Tlie CDHS data {E = 70 GeV), tliougfi not entering tfie fit, are also shown. 
The full and dotted lines show the f itl results. The dotted lines describe the 
data rejected from f itl [x < 0.1). The nuclear correction for these particular 
data are taken from ref. 0]. The error bars represent the quadratic sum of 
statistical and uncorrelated systematic errors. The additive bin constants are 
c{xi) = {12, 11, 10, 9, 8, 7, 6, 5, 4, 3, 2, 1.5, 1, 0.5} 
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Figure 14: Results of f itl for the longitudinal cross sections and structure func- 
tions. Comparison of R and Fl calculated in f itl with: a) the BCDMS and 
NMC measurements oi R — aL/cTT] b) the HI extraction of Fl] c) the CDHSW 
measurements (the dashed hne is F^, the sohd line is F^). None of these measure- 
ments enter our fits and only those corresponding to > 3.5GeV^ are shown. 
The error bars represent the quadratic sum of the statistical and systematic er- 
rors. Each point of these plots corresponds to different values of Q^. The curves 
are obtained by interpolating smoothly between the calculations performed at 
these points. 
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Figure 15: The strange distribution function of fitl, with its error band, for 
three different values of . The CCFR dimuon determination (full circles with 
statistical and systematic errors added in quadrature) is also shown. 
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Figure 16: Mean value of the pull distribution as a function of x (see text) for 
CDHSW anti-neutrino (a and c) and neutrino (b and d) data. Three fits are 
compared (see text). The arrows describe the changes of the pull values when 
going from one fit to another. 
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Figure 17: Results of f it2 for: a) the difference x{s — s) and b) the ratio s/s at 
g2 = 20GeV2. In the box b) the result of CCFR is also shown. 
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Figure 18: Difference between the uFe and uFe CDHSW differential cross- 
sections (see eq. (p6D) at fixed x and ~ 20 GeV^ as a function of YL = 
1 — (1 — Tlie solid line corresponds to f it2, the dashed line to f itl. The 
exact values of are 22, 21.6, 20.3, 18.6 GeV^ for x = 0.35, 0.45, 0.55, 0.65, 
respectively. 
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